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Abstract. We examine the scaling trends in particle multiplicity and flow observables between SPS, RHIC 
and LHC, and discuss their compatibility with popular theoretical models. We examine the way scaling trends 
between SPS and RHIC are broken at LHC energies, and suggest experimental measurements which can further 
clarify the situation 



1 Before the LHC: logarithms and triangles 

The azimuthal anisotropy of mean particle momentum ( 
parametrized by it's second Fourier component U2), thought 
of as originating from the azimuthal anisotropy in collec- 
tive flow ("elliptic flow"), has long been regarded as an im- 
portant observable in heavy ion collisions. The main rea- 
sons for this is that elliptic flow has long been thought to 
be "self-quenching" [ 1.2 1: The azimuthal pressure gradient 
extinguishes itself soon after the start of the hydrodynamic 
evolution, so the final vi is insensitive to later stages of 
the fireball evolution and therefore allows us to probe the 
hottest, best thermalized, and possibly deconfined phase. 

In addition, as has been shown in, elliptic flow is highly 
sensitive to viscosity. The presence of even a small but 
non-negligible viscosity, therefore, can in principle be de- 
tected by a careful analysis of 02 data. 

Indeed, one of the most widely cited (in both the aca- 
demic and popular press) news coming out of the heavy ion 
community concerns the discovery, at the relativistic heavy 
ion collider ( RHIC ), of a "perfect fluid", also sometimes 
referred to as "sQGP" (strongly coupled Quark Gluon Plasma) 
H3I4I5I6I7I . The evidence for this claim comes from the 
successful modeling of RHIC V2 by boost-invariant hydro- 
dynamics. 

Going further in our understanding is hampered by the 
large number of "free" (or, to be more exact, poorly under- 
stood from first principles) parameters within the hydro- 
dynamic model: While the equation of state is thought to 
be understood from lattice simulations 1 8 1, the behavior of 
shear and bulk viscosity is quantitatively not known around 
T c , where it is expected the temperature dependence could 
be non-trivial 1191 1 01 1 111 21 . The same goes for the large 
number of second-order transport coefficients. While we 
have some understanding of the initial transverse density 
of the system (its dependence on the transverse coordi- 
nate is thought to be either a "Glauber" superposition of 
p-p collisions 0141151 or a partonic semi-classical "color 
glass" |[T6l ). we do not as yet have control over the degree 
of transparency of the system, the amount of transverse 



flow created before thermalization (thought to be necessary 
to make the data agree with particle interferometry mea- 
surements [ 17 1), or of the interplay between the "medium" 
and the surrounding "corona" of peripheral p-p collisions 
H18I191 . A model incorporating "all physics", therefore, is 
expected to have a lot of correlated parameters which will 
be highly non-trivial to disentangle. 

A tool with the potential of overcoming these difficul- 
ties is scaling naturalness. Experiments have collected an 
extraordinary amount of flow data, encompassing a wide 
range of Energy ( Vs),centrality (parametrized by number 
of participants Npan), system size (mass number A of the 
nuclei), rapidity ifl particle species and transverse momen- 
tum (pj). 

The experimental data collected allows us to "scan" ob- 
servables dependence on variables relevant to the theory, 
and to see if the observable change when the same variable 
is obtained in different ways (for example, flow at mid- 
rapidity of a lower -\fs collision compared with flow at the 
fragmentation region of a higher -\fs collision,at the same 
multiplicity density dN/dy). 

This work qualitatively examines the scaling of multi- 
plicity with -\fs, and the scaling of flow observables with 
both yfs and dN/dy, between \fs = 19.6 GeV (SPS and 
lower RHIC energies) and V* = 2760 GeV (LHC ener- 
gies). 

This exercise, even when done in a very qualitative 
level, can yield significant insights into the dynamics since 
multiplicity and flow observables's scaling with rapidity 
and were remarkably simple [20 211221231 . The main 
trends examined in this work are: 

- Limiting fragmentation of dN/dy: At high rapidities, 
the slope of dN/dy becomes independent of reaction 
energy, the multiplicity curves, plotted w.r.t. y—ybeam 
if In y[s/2 can be superimposed 
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1 In this article we interchangeably use the rapidity y = 
tanh(p-/£) and the pseudorapidity rj = tanh(p-/p). Away from 
mid-rapidity and low p T the two are the same to a very good ap- 
proximation. Other than the label of Fig. 1, taken from 1 20 1, 77 in 
the text of this paper refers to viscosity, not to pseudorapidity 
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. Scalings of multiplicity and flow between SPS and top RHIC energies. Left panel shows limiting fragmentation, middle panel 
|arithmic dependence of multiplicity with energy, and right panel limiting fragmentation of elliptic flow. Experimental data taken 
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- Logarithmic dependence of multiplicity on -\fs: Mul- 
tiplicity at mid-rapidity scales approximately linearly 
with number of participants and logarithmically with 
the center of mass energy -\fs 

- Limiting fragmentation of V2'- the limiting fragmenta- 
tion can also be observed for elliptic flow. Consequently, 
central elliptic flow scaled by eccentricity also depends 
monotonically on dN/dy scaled by the area [23 1 

Fig. Q] summarizes these trends. 

The first of these phenomena was noticed a long time 
ago and is relatively straight-forward to understand [24|: 
In the co-moving frame of one of the nuclei, the other nu- 
cleus looks like a very thin Lorentz-contracted pancake. 
At high enough rapidity, the pancake is thin compared to 
all other scales of the system (basically the proton size 
Aq 1 cd ). Hence, dynamics should not be sensitive to how 
boosted the frame is w.r.t. pancake as long as the boost is 
"large". These concepts can be naturally implemented in 
Regge "string-based" models [25 26 27] and are compati- 
ble with the partonic description of hadrons [ 28 1 . 

The question arising naturally is... how large does "large" 
rapidity need to be? Looking at rapidity distributions, one 
can see that it can be very small indeed, even at ultra- 
relativistic energies: Limiting fragmentation breaks off into 
a rapidity plateau which, contrary to the predictions of Feyn- 
man scaling [29 30], is only Ay ~ 1 wide, a width that 
is largely independent of Vs. This fact provides a natural 
explanation for dN/dy ~ In V*: Since the total width of 
the rapidity distribution is ~ In by kinematics, and the 
width of the tip of the rapidity distribution is ~ ~\fs°, the 
height of the rapidity distribution should also ~ In V^- 

This apparently simple "explanation" hides a very non- 
trivial initial state dynamics: The reasoning used to explain 
limiting fragmentation, when extended to mid-rapidity, pre- 
supposes two "pancakes" passing through each other, each 
much thinner than Aq Cd . Hence, one would expect dynam- 
ics to be invariant under boosts, in other words a large (un- 
til y is a large fraction of yn m ) plateau around mid-rapidity. 

This is equivalent to the Feynman description of a boost- 
invariant multiplicity distribution |29 30) up to "high" ra- 
pidities Ay ~ O(yum) ~ 0(1) In s/s, a description natural 
in an asymptotically free theory, since an interaction mov- 
ing two high-rapidity partons into mid-rapidity has to be 
"hard" and hence suppressed. Instead, experimental data 



says that the rapidity pleateau is either non-existent or small 
(Ay ~ 1, of the order of the thermal smearing expectation) 
and independent of -\fs. The simplest parton-string models, 
where string ends predominantly originate from valence 
quarks, can not explain such a dependence, and hence can 
not describe dN/dy ~ In yfs 

The most natural way of reconciling limiting fragmen- 
tation away from mid-rapidity with the absence of boost 
invariance seems to be Il22l to maintain a "stringy" ini- 
tial state at formation time To, but admit that the space- 
time locus of collision (the spacetime region at zero rapid- 
ity),rather then the target and the projectile, act as sources 
of most string ends. A quantitative model of this type within 
QCD is, however, to date lacking, through something sim- 
ilar can be achieved by allowing for low x sea partons to 
form "small chains" (separated by very little rapidity) |25 1, 
or by representing these small —x partons as string excita- 
tions [26). 

If this scenario is correct, then, initially, flow is boost- 
invariant [29 30] (the spacetime rapidity is equal to the 
flow rapidity), but the initial density strongly depends on 
rapidity, and hence boost-invariance is not a good approx- 
imation, except very close to mid-rapidity. A generic pre- 
diction of this scenario is the breakdown of limiting frag- 
mentation when dN/dy stops growing logarithmically with 
In -\fs, something that can be studied at the LHC,as the next 
section will show BP . 

If, as commonly thought, elliptic flow is hydrodynamic 
in origin, the density rapidity variation means that at high 
rapidity most of the system will be hadronic, while at more 
central rapidity it will be partonic. Assuming the start of 
the expansion proper time T eq is only weakly dependent 
on -\fs (assuming otherwise generally breaks limiting frag- 
mentation [22 1), this seems to happen at y - yb eam ~ 1 Il22ll 
for asymptotically high energies. This region has been ac- 
cessed at both SPS and RHIC energies, and is shown in the 
distributions of Fig Q] 

Given these considerations, limiting fragmentation of 
V2 is particularly surprising: Assuming the rapidity depen- 
dence of flow observables is "encoded" in the initial den- 
sity (and associated intensive properties: T, rj/s etc.) rather 
than in the transverse size (S ~ A 2li at all y), and assum- 
ing subsequent evolution is local in y, we should expect 
the scaling of i>2 to follow the scaling of dN/dy provided 
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the intensive properties are either invariant with the parton 
density at equilibriump(Y e? ) (r eq is the time hydrodynamic 
evolution starts, ~ the mean free path \ m j p ~ tj/{Ts)), or 
change monotonically with p(j eq ), throughout the rapidity 
range. However, in a cross-over from a QGP to a hadron 
gas, intensive properties of the system should not change 
monotonically with p(j eq ). 

In Fig. [2] we summarize the expected changes: The vis- 
cosity to entropy ratio rj/s is expected to jump from the 
relatively high value of the hadron gas (77/s ~ N% in a gas 
of mesons and glueballs, where N c is the number of colors 
[32|) to the low value of strongly interacting QGP (77/5 ~ 
O (0. 1 ) [ 1 3 1), and then slowly increase to the asymptot- 
ically free weakly coupled QGP (77/ s ~ O (^~ 2 ) N®, where 
A is the 't Hooft coupling constant lfT3l ). In addition, the 
speed of sound is expected to have a dip in the cross-over 
region, whose depth is at the moment not well determined 

12m 

In a wide variety of models, elliptic flow 02 depends 
on the eccentricity e, and should be approximately propor- 
tional to it. This essentially follows from Taylor-expanding 
the solution of whatever dynamical equation obeys in e, 
since e is small and dimensionless, and since by symmetry 
with no e there is no elliptic flow. 

We also know that ^/e decreases if viscosity is turned 
on, i.e. if K increases. Hence, it is quite natural that, as 
suggested in B34I231 
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where K is the Knudsen number and Kq ~ 0(1) is a pa- 
rameter specific to the theory. Furthermore, the transverse 
Knudsen number at a given mean free path l m f p is 



K~ 



dN 



l m fpS dy 



(2) 



This formula assumes just boost-invariant flow, as well as 
a a time-scale [ 34] t„2 = VS/c s for the building up of v%, 
where VS (~ N l J a 3 rl , Not to be confused with the center 
of mass energy V* or the entropy density s) is the initial 
transverse size of the system (which, as we argued earlier, 
is independent of rapidity). It should be noted that going 
beyond this rough approximation for r„2 worsens scaling, 
since fl] t„2 rises and saturates with increasing density 
1 35 1 due to the self-quenching of elliptic flow. 

The derivation of Eqs.[T]and|2]follows straight-forwardly 
Il34l from density formula [ 30 1 



(3) 



1 dN 
P ST eq dy 

and Taylor expanding 

We further remember that ^| depends on the equa- 
tion of state, i.e. on the speed of sound. By a leading order 
expansion argument, and remembering that the asymptotic 
expansion speed of a Godunov-type hydrodynamic shock 
wave ~ c s 0611 . it can be seen that ^ ~ c s (numeri- 
cal simulations lend credence to this scaling, see |34|). We 
also remember that the mean free path l m f p ~ j-. 
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Fig. 3. (color online) The v 2 dependence on rapidity given initial 
conditions reproducing limiting fragmentation, and subsequent 
hydrodynamic evolution. The superimposed dashed line shows 
the pre-equilibrium (r = Td yn ) partonic density motivated by the 
dN/dy ~ In ^ and well reproduced by models such as [ 25 1 



Putting everything together, and neglecting the differ- 
ence between y and the pseudo-rapidity 77 (small in the 
fragmentation region away from mid-rapidity), we get that 



— ~ C s (T eq ) \ 1 



^(Ar^firr 1 ) 



Ts 



(4) 



we believe that when T > T c rj/s «: l,c s ^ 1/ V3, when 
T ~ T c c s <K 1/ y3 and r//s is at a minimum, and when 
T < T c c s goes back to a value not too different from 1 / y3 
but rj/s increases to > 1. (Fig. |2j. 

Additionally, close to mid-rapidity the plasma lifetime 
should » t„2, so V2 saturates and becomes independent of 
initial density (Eq.|4] over-predicts v^je). In the less dense 
region, however, the plasma lifetime < r„2, so V2 should be 
approximately proportional to the initial density (Eq.|4]is a 
good approximation). 

On the other hand, N par , should be independent of ra- 
pidity and pseudorapidity while the initial T(T eq ) should 
smoothly change as ~ (dN/dy) 1/3 OOll . 

We immediately see that the scaling seen in the right 
panel of Fig.Q]is not compatible with a modified BGK ini- 
tial condition, or indeed any initial condition without an 
unphysically finely tuned correlation between the size of 
the system and intensive parameters Il23l . In the suppos- 
edly long-lived ideal fluid mid-rapidity region, t^/e should 
be considerably flatter than dN/dy due to the self-quenching 
of i>2- At the critical rapidity where T eq ~ T c , vije should 
dip due to the dip in the speed of sound, and in the frag- 
mentation regions where Tq < T eq V2/e should go down 
more rapidly than dN/dy due to the rise in rj/s. The ra- 
pidity at which V2 vanishes should in general be different 
from the rapidity at which dN/dy does, due to V2 additional 
dependence on 77/ ^ and system lifetime. These expected 
trends are summarized in Fig. [3] 

These considerations have prompted us Il22l to doubt 
the hydrodynamic paradigm, and to search for a model 
where V2 ~ dN/dy independently of the thermodynamic 
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Fig. 2. Breakdown of scaling with temperature of the parameters of the hydrodynamic model 



properties of the system: Note that at formation time To 
,when dynamics starts, the system is partonic throughout 
the whole rapidity range So a far-from equilibrium system 
expanding very early will be immune from these consider- 
ations. However, hydrodynamics dictates that flow starts at 
local equilibrium T eq much later then the formation time. 
Only at T eq the system "knows" whether it is in the low 
rj/s partonic phase or in the hadronic phase. Perhaps the 
most straight-forward way of forming flow before equi- 
librium is to assume a Knudsen number ~ 1, which, as 
has been shown in [37|, can generate a significant amount 
of V2 (through not enough to describe experimental data): 
Since initial conditions are partonic throughout, therefore, 
V2 can be insensitive to weather the "equilibrium temper- 
ature" at that rapidity is above or below T c , thus ensuring 
that U2 depends monotonically on dN/dy [22. 38 1. Strong 
mean fields ~ dN/dy could then conceivably bring 02 up 
to experimental values without disturbing the scaling [39|. 

The trends described in this section allowed to predict 
both multiplicity and flow observables at the LHC 11201211 : 
It was hoped that this scaling would allow for a rigorous 
test of hydrodynamics, since the projected increase in flow 
would place vi above the hydrodynamic limit. As usual, 
however, data confounded facile predictions. 



2 Post LHC: More Multiplicity and same 
flow 

Experimental data has a way of spoiling simple and elegant 
trends! BP has shown that multiplicity has grown con- 
siderably faster than logarithmically, as ~ s , though, it 
should be noted, far slower than Landau hydrodynamics~ 
s 1/4 EU) and AdS/CFT Shockwave predictions ~ s 1/3 PT1 . 
suggesting that initial transparency is still very high, in line 
with the stringy models described in the previous section. 

V2, on the other hand, has grown approximately loga- 
rithmically 11421 : If one renormalizes the 02 triangle shown 
in |f7l I with the ratio of eccentricities in the relevant cen- 
trality bins (The eccentricity of 20-30% centrality event 
class of Pb-Pb reported in |g2) is ~ 20 - 30% higher than 
the eccentricity of 0-40% used in [20], as calculated from 
11431 and fl3]), it is clear V2 is in the region predicted by 
0201211 or slightly below it. 



Thus, the 1)2(77) scaling is probably in line with the trends 
from 02OI21I . but the multiplicity scaling of [20 1 is abun- 
dantly broken. This makes it likely that the integrated V2/e 
vs (\/S)(dN/dy) scaling, holding from AGS to LHC, is 
also broken. 

Looking at the preliminary HBT data 11441 . one con- 
firms the U2 results: The nice scaling of R ut,side,long with 
dN/dy 1 ^ 3 [45 1 is most likely broken due to the growth of 
dN/dy. While this is not shown explicitly, one can see [44| 
that RouiRsideRiong ~ In V* (as expected from the old scal- 
ing assuming dN/dy ~ In -\fs). however, dN/dy grows 
faster then In VA 

Its as if "flow" observables {u2 and HBT) maintain their 
old In V* scaling, but multiplicity grows faster. In this con- 
text, one can suggest two "extreme" scenarios: 

(i) The "soft" particle production still follows the scal- 
ings suggested in 1201211 . but the unthermalized mini- 
jets contribution to multiplicity can not anymore be ne- 
glected when predicting multiplicity. These minijets, 
however, are not in thermal equilibrium, and only "soft" 
particle dynamics dictates flow. 

(ii) The "extra contribution" above Il20l is fully thermal- 
ized and part of the system's soft dynamics. Flow scales 
with the total dN/dy, but it's scaling is more complex 
than that suggested by the right panel of Fig.Q] 

These, of course, are extremes, and reality can be some- 
where in the middle 0. A third possibility, that the extra 
multiplicity reflects the entropy created by the viscous evo- 
lution of the not-so-perfect fluid, seems unlikely in light of 
the good fit obtained of dN/dy as a function of centrality 
with the most popular initial state models [46 47 48 1; The 
overwhelming part of the entropy of the system seems to 
be there from the time of start of the dynamics To- 

As fig.|4]shows, a measurement away from mid rapid- 
ity, between the rapidity corresponding to SPS RHIC en- 
ergy, (respectively ~ 2.5 and ~ 1 .3) at the LHC is essential. 
As argued in the previous section, the simplest explana- 
tion for the In -\fs scaling of both multiplicity and V2 is that 
limiting fragmentation holds up to close to mid-rapidity. 
If this is true and scenario (i) holds, then, since integrated 

2 For example, 1271 seems to find that violations of limiting 
fragmentation happen at the same time for dN/dy and flow ob- 
servables 



Hot and Cold Baryonic Matter - HCBM 2010 



Scenario (i) 



dN/dy 



Scenario (ii) 

A 



; LHC 
; (inflection 
/ above RHIC 
' rapidity) 



dN/dy 



/ LHC 



RHIC 



y-y 

J -foe 



RHIC 



v2 
e 



Deam 
LHC 



J ,--fo( 



earn 



v2 
e 



/ RHIC 



-'LHC 



RHIC 



y-y 



earn 



y-y 

J -foe 



earn 



Fig. 4. The two interpretions of the trends between RHIC and LHC, and the expected consequences for the rapidity distribution in each 
scenario. Scenario (i) presupposes a "soft" background obeying the scaling of |20|, and an unthermalized extra "jetty" contribution. 
Because the jetty contribution is unthermalized, flow maintains limiting fragmentation (bottom panel). Multiplicity maintains limiting 
fragmentation up to rapidity smaller than that corresponding to ^ = 200 GeV (y = 1.3 ). In the other limit, limiting fragmentation 
is broken throughout the fireball's rapidity range due to the contribution of the "extra" multiplicity to flow |49 |. Reality could be, and 
probably is, in between these limits 



V2 is dominated by soft particles, limiting fragmentation 
of U2 should continue (not impossible, given the approx- 
imately logarithmic dependence of V2 in 11421 ). and limit- 
ing fragmentation of dN/dy should gradually break above 
the rapidity corresponding to RHIC energy. On the other 
hand, scenario (ii) would naturally predict a breaking of 
limiting fragmentation at all rapidities between RHIC and 
LHC, since the extra particles produced at mid-rapidity 
will generate longitudinal and transverse flow [49|. Exper- 
iment will tell us shortly which is the case. 

What do we make of the increase in U2? Contrary to 
earlier models, hydrodynamic codes with both viscosity 
and a long after-burner phase do allow for elliptic flow to 
rise above the (initially thought to be) "ideal fluid" limit 
of RHIC [50]. The compatibility of such a complex model 
with simple scalings, however, requires investigation. 

To start interpreting these results within the hydrody- 
namics paradigm, examining a further scaling is necessary: 
The differential scaling of t^Cpr)- As can be seen when 
comparing |51 1 with [42 52], below -\fs ~ 62 GeV, v%{pt) 
saturates, and the increase of integrated U2 is driven not 
by an increase of viipr) but by a higher average (pr)- An 
ideal hydrodynamics scan [£35 1 with different energies pro- 
duces just such a behavior, provided rj/s and T f reezeou , do 
not change across the considered V^- A similar simula- 



tion with viscous hydrodynamics [53|,however, shows that 
such scaling very quickly breaks down in case rj/s changes. 

An interpretation of these findings [52] is that rj/s sat- 
urates at - 62 GeV and stays constant up to LHC 
energies. But then, what does one make of limiting frag- 
mentation of V2 up to V* — 19.6 GeV, and of the scaling 
of u 2 /c vs (l/S)(dN/dy) scaling between AGS and RHIC 
energies? 

Alternatively, it can be postulated that changes in the 
initial temperature, rj/s and Tf reezeol „ somehow compen- 
sate each other in the final flow observables. Weather this 
is possible without fine-tuning is still an open question 
1 50 1 . It should be noted, however, that, while changes in 
Tfreezeout are reasonable in hydrodynamics (Since freeze- 
out happens when the mean free path l m f p ~ R, the system 
size, it is reasonable to assume that the higher the initial 
energy, the lower Tf reezeout ), the data we have seems to dis- 
allow such changes over the available -\fs: If this was the 
case, one would expect a breaking of the scaling of HBT 
radii with dN/dy (due to longer lifetime in higher density 
events), as well as a depletion of ratios of particles such as 
K* /K with increasing density 0541551 . It is fair to say no 
such systematics exists [45 56 57 1. 

If, as proposed in 11521 . the V2(pr) scaling is the "fun- 
damental" one, the critical observable to look for here is 
(Pt) ( V*) : Does (pr) (thought, unlike U2 to increase uni- 
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formly throughout the lifetime of the system) scale, as i>2, 
with In -\fs, or does it scale with dN/dy (faster than In -\fs 
at LHC energies)? Only in the latter case can we be sure 
that the increase of total 02 is dominated by (pr) changes. 
The latter case is also favored in scenario (ii) (Scenario (i) 
would still prefer a parametrically slower growth for (pr) 
wrt dN/dy if out-of equilibrium pj > 2 - 3 GeV parti- 
cles were cut out), and natural in saturation-based scenar- 
ios l58t 

Independently of these considerations, the problems re- 
lated to describing V2U)) with hydrodynamics will remain. 
If, at moderate (corresponding to RHIC V?) rapidities lim- 
iting fragmentation of V2 is restored, looking for {pj) and 
HBT radii dependence of rapidity will clarify to what ex- 
tent are flow properties correlated with the equation of state. 

If flow properties (u2, (pr) and HBT radii) continue ex- 
hibit limiting fragmentation at moderate rapidity, models 
of flow generation far from equilibrium such as M37I39H 
will need to be given serious consideration. 

In conclusion, we have reported, and attempted to in- 
terpret, experimental scaling trends of multiplicity and flow 
properties in heavy ion collisions at SPS, RHIC and LHC 
energies. Between SPS and RHIC energies both multiplic- 
ity and flow observables ~ In -\fs and exhibit limiting frag- 
mentation in rapidity. While the multiplicity scaling can 
be naturally described by popular string-based approaches, 
flow scaling presents a challenge to the generally accepted 
hydrodynamic model. 

Intriguingly, at the LHC, while flow observables seem 
to scale with In yfs, multiplicity grows faster. This could 
mean that a "hard" (possibly-nonflowing) component to 
the "soft" (and flowing) background becomes non-negligible. 
Determining whether this "hard" component flows or not 
is therefore important. Examining limiting fragmentation 
of i>2 , (pr) and HBT radii, and seeing how {pj) grows 
with i{s, can shed light in this direction. 
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